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^ '. Abstract 

We derive a QCD sum rule for the flavour-singlet axial coupling constant g A from a two point 
correlation function of flavour-singlet axial vector currents in a one-nucleon state. In evaluating 

<n : 

the correlation function by an operator product expansion we take into account the terms up to 



dimension 6. This correlation function receives an additional two-loop diagram which comes from 
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an (anti-)instanton. If we do not include it, g A is estimated to be 0.8. However, the additional 
diagram due to instantons contributes negatively and reduces g A towards the experimental value. 
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Axial coupling constants are defined by nucleon matrix elements of axial currents at 
zero momentum transfer. Since an axial current, q{x)j IJ ,'j5q(x) , is a spin operator, the 
flavour-singlet axial coupling constant g A represents the fraction of the nucleon spin carried 
by quarks. In the naive parton model, g A is expected to be close to 1. However, an 
unexpected small value of g A was found from the EMC experiment, which implies the 
quarks contribute only a small fraction to the proton's spin. This has led to the so-called 
"spin crisis" and raised a number of questions of understanding the dynamics of the proton 
spin P . A number of subsequent experiments have been performed. The results are in the 
range g A = (0.28 — 0.41), see [2] for a recent review. 

The investigations of g A by QCD sum rules have been done so far by the authors in 
Refs. |8l9j . Ioffe and Oganesian |H] derived a QCD sum rule for g A by considering a two-point 
correlation function of nucleon interpolating fields in an external flavour-singlet axial- vector 
field. In their sum rule, g A is expressed by expectation values of QCD composite operators 
induced by the external field. An important contribution comes from dimension 3 term 
in the operator product expansion, which is related to the derivative of QCD topological 
susceptibility x'(0)- They found the lower limit of g A and x'(0) from the self consistency 
of the sum rule: gf > 0.05, x'(0) > 1.6 x lO^GeV 2 . 

The authors in Ref.jU] considered a three-point function of nucleon interpolating fields 
and the divergence of a flavour- singlet axial- vector current. They took into account chiral 
anomaly by using the anomaly relation from the very beginning. The form factor, g A (q 2 ), is 
related to the vacuum condensates of the quark-gluon composite operators through a double 
dispersion relation. To know g A (q 2 ) at q 2 = one must evaluate the correlation function 
at zero momentum. Although the method to evaluate it is known [10] . it involves large 
uncertainty. 

Recently, we have proposed a new method to construct QCD sum rules for axial coupling 
constants from two-point correlation functions of axial- vector currents in a one-nucleon state 
[TT] . With the method, the axial coupling constants are expressed in terms of the it- N and 
K-N sigma-terms and the moments of parton distributions. We have seen good agreement 
with experiment for the non-singlet constants, g A and g A . 

In this paper, we extend the previous work to the case for g A . For the calculation of 
g A we need to fully take into account the chiral anomaly. Since the origin of the chiral 
anomaly is considered to be instantons, one might suspect the anomalous suppression of g A 
is somehow related to instantons. We therefore evaluate their effects on g A . 

Following Ref . ^^ 5 we consider a correlation function of flavour- singlet axial- vector cur- 



rents in a one-nucleon state: 

nS(g;P) =iJ*xe*(T\j®(x)J%\0))) N , (1) 

where g M = (u>, q) and the nucleon matrix element is defined by (. . .}jv = 
(V 2 ) Es [(N(PS)I . . . |iV(P5)) - (. . .)o(JV(P5)|iV(P5'))], where P» = (E,P) is the nucleon 
momentum (P 2 = M 2 , M is the mass), S the nucleon spin, (. . .) = (0| . . . |0), and the one- 
nucleon state is normalized as (N(PS)\N(P'S')} = (2tt) 3 (5 3 (P -P')5 S s'- The flavour- singlet 
axial- vector current is defined as 

3% ( x ) = u(x)^^u(x) + d(x)^^d(x) + s(x)7 At 7 5 s(x), (2) 

where u, d and s are the up, down and strange quark fields, respectively. 
We write a Lehmann representation of Eq.(jT|): 

rr(0)/ td\ f°° a //VO^''^-^) (r>\ 

njJ(w,q;P)= / du — , (3) 

J~oo UJ — UJ 

where p^ v is the spectral function. We derive a Borel sum rule from Eq.(j3J) for the even 
function part of Eq.tJTJ) in uj, Ulu! (^ 2 , q; P)even, as [TT] 

/oo 
du'u'exp (-lu' 2 /s) p$(u', q; P), (4) 

-co 

where s denotes the square of the Borel mass and B the Borel transformation with respect 
to uj 2 . In Eq.pjl the left hand side is evaluated theoretically, which give rise to a Borel 
transformed QCD sum rule. 

Let us now consider the physical content of the spectral function with the insertion of 
intermediate states between the currents. Here the lowest one is a one-nucleon state. The 
continuum state consists of ^'-nucleon states, excited nucleon states and so on. There is 
an energy gap between the nucleon pole and the continuum threshold. The contribution of 
the one-nucleon state to the spectral function is expressed in terms of axial coupling con- 
stants, because the nucleon matrix element of j^ 5 is written as (N(PS)\j fl5 (0)\N(P'S')) = 

u(PS) 5 , V(? 2 )7/i75 + ^A / (? 2 )?//75 u(P'S'), where u(PS) is a Dirac spinor and q = P' — P 
jllj . The contribution of the continuum state becomes small in the Borel sum rule, since 
it is exponentially suppressed compared to that of the one-nucleon state because of the en- 
ergy gap. Therefore, it is allowed to use a rough model of the continuum: The form of the 
continuum is approximated by the step function with the coefficient being the imaginary 
part of the asymptotic form of the correlation function in the OPE ^H]- In the present 
case, however, the continuum contribution to the spectral function is absent within the ap- 
proximation, because the perturbative part is subtracted from the definition of Eq.(fTj). This 



means that the continuum contribution may be very small at least in the high energy region. 
We therefore neglect the continuum contribution in this work. 

Hereafter we consider the correlation function in the rest frame of the initial and final 
nucleon states and contract the Lorentz indices of the currents. Expanding the right hand 
side of Eq.(j3J) in powers of |q| 2 , we find the coefficient of |q| 2 is proportional to |gV(0)| 2 
fTT] , h A (0) contributes to higher order terms since h A (q 2 ) has no singularity at q 2 = 0. 
From the first derivative of Eq.(@J) with respect to |q| 2 we obtain the desired QCD sum rule 
at |q| 2 = : 



aB[n(°V, q ) cvcn ] 



<9|q| 



|q| 2 =0 



4ia 2 - (5> 



where n(w, q) = 11/ (w, q; M, 0). 

Let us now turn to the evaluation of IT( )(g). n(°)(g) consists of the following two terms: 

q=u,d,s q,q'=u,d,s 

where C(q) q is a connected or "one-loop" term which is given by 

C(q) q = -iJ d\x - y)e^^ (Tr {T [q(y)q(x)\ ^^T [q(x)q(y)\ ^ l5 }) N (7) 

and D(q) qq i a disconnected or "two-loop" term: 

D(q) qq , = i [ d\x - y)<&*-T) (Tr { lfll5 T [q(x)q(x)}}Tr {f^T [q > (y)q'(y)}}} N . (8) 



Note here that only the correlator of flavour singlet currents receives the contributions of 
the "two-loop" terms. Indeed, for the correlation function of iso- vector currents, j § = 
(l/2)(u , j f j j '-f 5 u — d'j f j i '-f^d), "two-loop" terms cancel with each other if we neglect the differences 
of quark masses: 



n (3) (<?) = iJd A xe^(T\j^(x)jS(0)]} N 



(1/2)" [C u (q) + C d (q) + D(q) uu - 2D(q) ud + D(q) 



dd\ 



= (l/2Y[C u (q) + Cd(q)}. (9) 

Similarly, for the 8-th component of flavour-octet currents, j^ 5 = (l/2\/3)(u , y ll 'y 5 u+d'y ll , y 5 d— 
2s7 M 7 5 s), 

n (8) (<7) = ijd 4 xe^(T\j^(x)J^(0)]) N 

= (l/2V3)' l [C u {q)+C d (q)+4C s (q) 

+D{q) uu + D{q) dd + AD{q) ss + 2D{q) ud - AD{q) us - 4D(q) ds ] 
= (l/2Vt) 2 [C u {q)+C d {q)+AC s {q)}. (10) 



We evaluate Eqs.(J7J) and ((HI) by a standard operator product expansion (OPE). Let us first 
consider the "one-loop" terms. In the OPE, operators of the leading terms are of dimension 
4. We take into account the terms up to dimension 6. The result is in the following: 



C{q) q = -^m q (qq) N - — (-^G 2 ) —i (qS-f IM D u q) N 
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—m q (qq) N - — I —G ) 

q 2 2q l \ it I n q 4 

22na 
-^ (q^X a q (u^u + d 1 ,X a d + s 1 ,X a s)) N 

p-!_ (S (q^X a q) (u lv \ a u + d lu \ a d + S lu X a s)) N 

32q»q v q X <f 

-t{qS>y fl D I/ DxD a q} N (11) 



qs 

where D^s are covariant derivatives, G 2 = G c \ v G a ^ v , and S denotes a symbol which makes 
the operators symmetric and traceless with respect to the Lorentz indices. 

We now discuss about the nucleon matrix elements in Eq.ljllJI. It is known well that 
m q (qq)N is related to the tt-N or K-N sigma-term as (m u + md)({uu)js[ + (dd)x) = 2Y*nN 
and (m s + m u )((ss)N + (uu)n) = 2T,kn- ((o: s /7i)G 2 )n is expressed by the nucleon mass and 
m q {qq) N through the QCD trace anomaly:((a s /7r)G 2 )Ar = -(8/9) [M - Y, q =u,d,s m 9^ N )- 
The matrix elements which contain covariant derivatives are related to the parton distri- 
butions as (Sqj lil D IMl ---Dn n q((j, 2 )) N = (-i) n ~ l A q n (ii 2 )T^^..^ n , where A n (fi 2 ) is the ra-th 
moment of the parton distributions at scale /i 2 , and T^^...^ = S [-P Ml -P M2 • • ■ P^ n ]- For the 
matrix elements of four quark operators, we apply the factorization hypothesis ^21 : the 
matrix elements are factorized by assuming that the contribution from one nucleon state 
dominates in the intermediate states: (0\02)n ~ {0\)n(02)q + {Oi)q{02)n- We apply this 
hypothesis to the following type of the nucleon matrix elements, which appear in Eq. ffTTj) : 
(qn^ a qfqrlu\ a qf)N = -(8/9)g^(qfqf)o(qfqf)N$f,f>, where / and /' are flavor indices. 

Let us next consider the "two-loop term", Eq.fjHJ). There is no contribution from this term 
as long as we do not account for higer dimensional terms (larger than dimension 6) in the 
OPE and perturbative corrections. However, "two-loop term" receives a contribution arising 
from instantons. We evaluate it under the "dilute instanton gas approximation" (DIGA) 
|17j . Namely we assume there exists only one instanton or anti- instanton in vacuum. We 
have two reasons why this approximation is expected to be valid. The first is that since we 
use the framework of QCD sum rules we are interested only in the short distance behavior of 
the correlation function. The second is that instantons in vacuum is sufficiently dilute [T7j . 
Indeed, it is known that the density of instantons is about 1 fm -4 and the critical size of an 
(anti-)instanton p c is p c ~ 0.3 fm. The value of p c is significantly smaller than the typical 
separation between instantons. 

In order to evaluate Eq.(|8~|). we first consider the correlation function in nuclear matter 
with its baryon number density being p#. Then Eq.fjHJ) is obtained as the first derivative 



of the correlator in nuclear matter with respect to ps, because, in general, expanding a 
vacuum expectation value of an operator O at finite baryon number density in powers of 
Pb the coefficient of the linear term is nothing but the nucleon matrix element: (0) PB = 
(O)o + Pb{0)n + ■ ■ ■■ Thus Eq.(jHJ) can be written as 
d 



D(Q) 



qq' 



d \ x _ y)jQ(*-v) (Tr { 7m75 T [q(x)q(x)]}Tr { T " 75 T [q\y)q\y)]}) 



PB 



D(Q) qq > = J dp 



MPiPb) 



PB=0 



l d PBj J PS= 

(12) 
where Q is an Euclidean momentum defined by Q 2 = —q 2 and (■ ■ ■) means that averaging is 
performed over all gauge configurations with the weight function exp(-S), where S is the 
Euclidean action. Under DIG A Eq. ()12|) becomes |2*T] 

~_d_ 

pp B 

d A xe lQx Tr [j^ 5 S q (x, x)] J d'ye'^Tr [ Jvl5 S q ,(y, y)\ , (13) 

where S q (x,y) is a quark (q) propagator in the field of an (anti-)instanton. We performed 
integrations over the center of an (anti-)instanton z and instanton size p with a weight 
ti(p,pb)- n(p, pb) is the tunneling rate at finite baryon number density. 
S q (x,y) is expressed as follows: 

^) = £^, («) 

A 

where ip\(x) is an eigen function of Dirac operator with the eigen value A: lfiip\(x) = Xip\(x). 
Then the dominant contribution comes from the zero-mode, ipo(x). In Eq.([13|). however, 
quarks do not propagate in zero-mode states. The reason is very simple. A zero-mode 
changes its chirality in passing through an instanton (see FigO] (a)). On the other hand, 
quarks created by an axial current have same chirality (FigQ (b)). So zero- modes are 
not allowed in Eq. (jl3j) and only non- zero- modes contribute. Non-zero-mode propagator 
Sg M (x,y), in which all non-zero-modes are summed up, 

Sf"(«)-E^ (15) 

satisfies the equation 

{# + im q )S^ ZM {x, y) = *«(* - y) - M*)4(v)- ( 16 ) 

Subtracting the zero- mode contribution in the right hand side, the remaining Si WZM (x,y) is 
ensured to be orthogonal to the zero-mode [22]. The solution of this equation, in general, 
has the following form |2*2"] : 

S» ZM (x,y) = ($ x - imq )A(x,y)^^ + A(x,y) 5 y ^^, (17) 



where D^ and D^ are covariant derivatives: D^=d^ —i^A^ and -D M = — <9 M —i^A -^ with 
r a 's are Pauli matrices and A a being an (anti-)instanton solution : 

a;(x) = 2^^4. (18) 

' x 2 x 2 + p 2 

Eq.(|18|) is the solution in the singular gauge. Hereafter we work in this gauge. An anti- 
instanton solution is obtained by replacing f\ a ^ v to r\ apbV . r\ aiW and fj atiu are the 'tHooft symbols 
fT7] , A(x, y) is the propagator of a scalar field which satisfies the equation 

(-D 2 + m 2 q )A(x,y) = 6^\x-y). (19) 

It is known that this equation is solved by 

A '^)-4^r^( 1 + fl" /2 ( 1 + S" /2 ( 1 + p2ll #)- (20) 

where r^ = (r,=pi) |22I23| . Here we have neglected 0{m 2 ) terms. The propagator in the 
field of an anti-instanton is obtained by interchanging r + and r~. 

The calculation of Tr \j^ b S q ZM {x^x)\ in Eq. (|T3"j) needs some care, since S^ ZM (x,y) is 
not well defined at x — > y. This is because S^ ZM (x, y) at x ^ y is not gauge invariant. The 
trace should be defined as limit of an gauge invariant expression in which a path ordered 
product is inserted: 



rx+e/2 T a 

* / —A a v {z)dz v 

Jx-e/2 A 



Tr [ lf , l5 S» ZM (x, x)} = lim Tr I lf , l5 S» ZM (x - e/2, x + e/2)exp 

(21) 
where the right hand side must be averaged over all the direction of the four vector e M . As 
a result we obtain 

The Fourier transform of this equation is given by 

J d 4 xe^ x Tr [j^ 5 S(x, x)] = -ip 2 Q^K (Qp). (23) 

Here i^o(^) is an 0-th modified Bessel function. Then Eq.([13|) reads 





D qq ,{Q) = j dp 



dp 



B 



Mp,pb) 



Q 2 p*K {Q P y. (24) 

Jps=0 



An important quantity in Eq. (|24|) is the tunneling rate, n(p, ps)- In normal vacuum, 
where vacuum condensation does not exist, the tunneling rate at one- loop was first given in 



Ref.|19j. After the pioneering work, it is now available in two-loop renormalization group 
invariant form |2(Jj . 

(m q p)(pn r^^r- (25) 

q=u,d,s 

where n (p) is that for quarkless theory: 



d^v ( 2n \ c ( 2n \ , , a „ </a <Ma .^HsW 



2N C 

n (p) = ^f ( -^— ) exp ( ^— ) {pp^^-^r-^ , (26) 

2e 5 / 6 
dm = ^(iV c -l)!(iV c -2)! eXp( ~ 1 - 511374iVc + °- 291746n /)> ( 27 ) 

/3o = yiVe - \n f , (28) 

A - T^ 2 - (t n - ~ K) »* (29) 

iV 2 -l 
70 = 3^-. (30) 

In Eqs.(|25| ) M|30 |) . yU is some arbitrary normalization point, N c = 3 and rif = 3 are the 
numbers of color and flavor. 

In physical vacuum, where vacuum condensation exists, according to Shifman et al. J21] , 
the current quark mass in Eq. (f25|) is substituted by a dynamical one. Shifman et al. con- 
sidered the tunneling rate for small size instantons. Then it can be expanded in powers of 
p. The coefficient of each term of the expansion should be a quantity characterizing the 
vacuum structure. They found that the the tunneling rate is written as a vacuum expecta- 
tion value of an "effective Lagrangian" , which has form analogous to a standard operator 
product expansion. The Lagrangian is as follows: 



ac = Mpx^r 70 ^ j n (™«p-^w?) 



q=u,d,s 
° • H 2„3 



H -it p 

32 \3 ' 



A ^ a - a X { 4 



uT a udT a d - ^uY a uv udT a uv d ) ( m s p - ^n 2 p 3 sTs 
9 4 



+ — • -7T p dabcuT^udT ^dsV 2 's + (2 permutations) 



40 3 

320 V 3 / 256 V 3 / ^ ^ 7 w 



(31) 



where T = (1 - 7s )/2, T a = (1 - 7s )/2 • (A a /2), T^ = a^(l - 75 )/2 ■ (A a /2) and / afec and 
<i a fe c are SU(3) symbols. For anti-instanton replace 1 — 75 to 1 +75. Then the tunneling rate 
at finite density is given by the average of Eq. (|31j) over the ground state of nuclear matter: 

n(p,p B ) = (AC) PB . (32) 



The expectation values of multi quark operators in Eq. (|52*|) are evaluated by applying fac- 
torization hypothesis. In this approximation, all the terms in Eq. (l31|) containing 75, a pu , 
X a /2 drop off [21]. As a result, Eq. (|32j) is reduced to the same form as Eq. (J25|) but with the 
current quark mass replaced by the "effective mass" : 



n{p,p B ) = n (p)(ppo) n/7 ° 
where the "effective mass" is defined by 



"ms^o) 



n m *q(P'PB)p, 
q=u,d,s 



(33) 



m q ~ -7T p (qq) PB . 



(34) 



In order to know the derivative of Eq. ()33|) with respect to p# in Eq. (J24JI , we must know the 
Pb dependence of (qq) PB in the effective quark mass. (qq) PB is expanded in powers of ps as 
(QQ)pb = (v<l)o + ^t~ d PB + ■■■ for q = u,d and (ss) PB = (ss) + y^t:pB + • • • [El, where 
y = 2(ss) n / ((uu) n + (dd)^) is the strangeness content of the nucleon. Then we obtain the 
derivative of Eq. (}3"3"j) with respect to p B as 



_d_ 
dp 



n(p,PB) 



B 



-'■kN 



PB=0 



n (p)(ppoT n °- 

3 m u + m d 

x K(p, 0)m* a (p, 0) + m* d {p, 0)m*(p, 0) + ym* u {p, 0)m*(p, 0)] 



(35) 

Now we have all the ingredients for deriving the QCD sum rule. Collecting all the terms 
in Eq.©, namely, the one- loop term Eq. (jllj) and the two- loop term Eq. (j24j) with Eq.([35|h 
and substituting Eq.© into Eq.©, we obtain the QCD sum rule for g A as follows: 



Wa 



(0)1 2 



M 

M 
+ — 
s 



j ttN 



28 

t 



1 - 



m s 



m u + m d 



+ 



■•KN 



56 



m. 



3 m s + m u 



- - 7 (A u 2 (p 2 ) + A d 2 (p 2 ) + A*(p 2 )) 



4:na s (qq) 



15M 3 



352 £ 



71 N 



27 m, 



Wd 



4ira s {ss) 



176 / 2S 



KJV 



->ttN 



[At(p 2 ) + Ai(p 2 ) + At(p 2 )]+I(s) 



where (qq)o = (uu)o = (dd) . In Eq.fj36| we assume m h 
contribution which is given by 



27 \m s + m u m u + m^ 

(36) 

rrid- lis) is the instanton 



Ks) 



9 / dp 
'0 



<9 



<9p 



B 



n{p,p B ) 



p' I dt(s — p cosh £ • s )exp(— p cosh £ • s). (37) 
J PB =o Jo 



Here the integration over p has been cut off at the critical size p c . We note that the value 
of p c estimated by Shuryak is close to the upper boundary for the validity of the expansion 
Eq.flHU) J2H- 



We show in FigQ the square of the Borel mass, s, dependence of \g A | in Eq. f|HU|) . In 
plotting the curve in Fig^ we used the following values of the constants in the right hand 
side in Eq. (j36)l . The n-N sigma-term is taken from Ref . |14j . which are E^ = 45MeV. The 
quark masses are taken to be m u = m^ = 7MeV, m s = HOMeV [T2] . Using the above 
values and the ratio y = 2(ss)n/((uu)n + (dd)^) = 0.2 given in Ref.[Hj, we can calculate 
the K-N sigma term averaged over the iso-spin states and the result is T>kn = 226 MeV. 
We calculated the moments of parton distributions adopting the LO scheme in Ref.|16j: 
^(lGeV 2 ) + ^(lGeV 2 ) = 1.1, ^(lGeV 2 ) + Aj(lGeV 2 ) = 0.13, A*(lGeV 2 ) = 0.03, 
^(lGeV 2 ) = 0.002. The vacuum condensates are taken from Ref. [12]. which are (qq)o = 
(—225 MeV) 3 and (ss)o = 0.8(qq)o. The normalization point /io in Eq. (j33j) was taken to 
be 1 GeV, which is the relevant scale for QCD sum rules. For the critical size of an (anti- 
)instanton we used the value p c = 0.3 fm. The upper curve in Fig |21 corresponds to g A 
in Eq. ([3l)j) but without the instanton contribution I(s). We see that the curve is well 
stabilized. In the stabilized region g A ~ 0.8. This is consistent with the well known fact 
that g A is about 30% suppressed due to relativistic effect compared with the naive quark 
model's expectation, but is much larger than the experimental value g A = (0.28 — 0.41). The 
lower three curves correspond to g A including the instanton contribution I(s) with different 
choices of p c . We see the Borel curve is extremely sensitive to p c and not stabilized. Therefore 
we cannot predict the value of g A . However, we can say that apparently instantons tend to 
lower g A compared with that without the instanton contribution. 

In summary, we have derived a QCD sum rule for g A from a two-point correlation 
function of flavour singlet axial-vector currents in one-nucleon state. In deriving the sum 
rule, we evaluated the correlation function by an OPE up to dimension 6. We have also 
took into account an additional contribution arising from an (anti-)instanton and evaluated 
it under DIGA. When we do not include the instanton contribution, g A is not so suppressed 
as the experimental value and is about 0.8. Including the instanton contribution g A tends to 
be suppressed compared with the result when we do not include the instanton contribution. 
Recently, Schafer and Zetocha J2H] have computed the axial coupling constants of the nucleon 
using numerical simulations of the instanton liquid. They found the isovector axial coupling 
constant is g A = 1.28, in good agreement with experiment, while flavour singlet coupling 
is g A = 0.77. g A comes from a connected part and OZI violating disconnected part of the 
three point correlation function. Taking into account only the connected part they found 
9a = 0.79, while the disconnected part is very small, g A (dis) = —(0.02 ± 0.02). It would 
be interesting if we can clarify the relation between the present result obtained by adding a 
single instanton contribution to the OPE-based QCD sum rules and that by the instanton 
liquid model. 

The author would like to thank Professor Y. Kondo and Professor S. Saito for valuable 
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FIG. 1: (a) Zero-modes propagating through an instanton. The solid lines correspond to zero- 
modes with the chirality left-handed (L) and right-handed (R). An instanton is shown as an open 
circle with "I". A zero-mode changes its chirality when passing through an instanton. (b) The 
diagram corresponding to the two-loop term, Eq.((8J). Quarks created by an axial current have 
same chirality. So quarks do not propagate in an zero-mode state, which changes its chirality in 
passing through an instanton. 



< 



t / 

i / 

i / ,' 

i / 

i / - 



0.8 

0.6 

0.4 

0.2 


0.8 1 1.2 1.4 1.6 1.8 2 

s (GeV 2 ) 

FIG. 2: The square of the Borel mass, s, dependence of \g A \ in Eg. (|36|) with and without the 
instanton contribution. The upper curve (long-dashed line) corresponds to that without the instan- 
ton contribution. The lower three curves show those with the instanton contribution for different 
choices of p c ; The solid line corresponds to p c = 0.3 fm, the long-dashed line to p c = 0.29 fm and 
the short-dashed line to p c = 0.31 fm. 



